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1. Introduction 

It is an interesting and long standing problem to understand whether the flavour-singlet axial 
U(l) symmetry is restored or not above the finite temperature transition in the chirally symmetric 
phase of QCD. At low temperature it is well known that the chiral symmetry is spontaneously 
broken while the axial U ( 1 ) symmetry is broken at quantum level by the presence of configurations 
with non-trivial topological structure. The semi-classical configurations that contribute to the axial 
charge are called instantons and it can be shown that zero modes of the Dirac operator are related 
to the presence of these and other configurations that have non-trivial topological charge. Hence, 
topology plays a major role in this subject. One clear example is the Witten-Veneziano relation [1, 
2] that connects the topological susceptibility of pure gauge theory to the heavy mass of the tj'(958) 
particle, the candidate to be the would-be Nambu-Goldstone boson of the axial U(l) symmetry. 

Our purpose is to study the behavior of meson correlators at high temperature, looking for 
signals of (approximate) degeneracy in all the singlet and triplet channels (a, 8, 7t,ri') toward the 
chiral limit. This would be an evidence of effective restoration of both chiral and axial symmetries, 
several intermediate steps are to be considered to control the systematic errors. The best avail- 
able discretization of the Dirac operator is the Overlap fermions [3]. Dynamical simulations with 
Overlap fermions are possible with current machines and were performed by JLQCD collabora- 
tion in the past years [4]. The price payed in order to be able to have a very precise fulfillment 
of the Ginsparg- Wilson relation [5] on the lattice is to fix the topology. Changing topology in a 
HMC trajectory requires generating rough configurations associated with near-zero modes of the 
hermitian Wilson operator (Hw)- A discontinuity arises in the HMC when one of these modes 
crosses zero. Treating this discontinuity is a costly task, although it could be achieved using the 
reflection/refraction methods [6]. One solution to this problem is to suppress near zero modes of 
Hw and to avoid the zero crossing, with the disadvantage of preventing topology change during the 
simulation. The simulation is no more ergodic, and some technique has to be developed to obtain 
the physical results at 6 = 0. This program was indeed developed by JLQCD collaboration (see 
[7] for details) and it works nicely at zero temperature. In order to accomplish our main target, we 
need to be sure that the same methods can be applied even in the finite temperature case. The first 
step towards the understanding and controlling the effect of fixing topology at finite temperature 
is to reproduce the known results of pure gauge theory. In particular we focus on the most rele- 
vant quantity, i.e. the topological susceptibility %,, and see if it can be reproduced with the same 
methods adopted at zero temperature. 

The rest of the paper consists of two main sections. The first is devoted to the discussion of 
our study of topological susceptibility in pure gauge theory at finite temperature. The second is the 
measurements of meson correlators, in flavor singlet and triplet channels, and the spectral density of 
the Overlap Dirac operator, where we show some evidences that the axial U ( 1 ) symmetry breaking 
is suppressed after the chiral phase transition. In the last section we will draw some conclusions. 

2. Pure gauge simulations 

In a previous paper [8], the JLQCD-TWQCD collaboration calculated the topological suscep- 
tibility at zero temperature in the case of two flavors of overlap fermions. They demonstrated that 



2 



Topological susceptibility and axial symmetry at finite temperature 



Guido Cossu 



it is possible to obtain physical quantities even with simulations at fixed topological sector. They 
checked the prediction of #PT, i.e. % t = (m q L)/Nf + 0(m q ), finding the expected linear behav- 
ior of %t versus the sea quark mass. The methods described in [7] were applied in that case, i.e. 
the topological susceptibility is measured by the long distance value of disconnected correlators 
of pseudoscalar operators (they are non-zero because in a fixed topology environment clustering 
property of field theory is violated). We will investigate numerically whether we can apply the 
same methodology at finite temperature too, with a preliminary test in pure gauge simulations, in 
order to be trained in dealing with this systematic error in the most interesting full QCD case. 

The correlators are calculated assuming that the major contribution comes solely from the 
lowest 50 eigenmodes. We tested for one value of j8 that this is really the case by changing the 
number of eigenmodes to 30 and 40. The result is that the long distance behavior saturates at 40-50 
eigenvalues, the short distance correlator has a worse approximation, as expected. 

In the pure gauge theory,there is an anomalous contribution to the pseudo-scalar meson corre- 
lator in the flavor singlet channel, which is called hairpin diagram [9, 10]: 

H(p)=f P ^- 1[ m 2 -z^rfp. (2.1) 

This is a zero temperature result, and we assume that it is valid at finite temperature with all param- 
eters depending on temperature. We actually find that this fits our data very well. Also, the results 
are independent of the valence quark mass, as it should be in the case of pure gauge theory. For 
statistical purposes, results are obtained with a joint fit of connected and disconnected part. The 
joint fit assumes an identical decay rate for the two correlators, long distance of connected part 
going to zero, to constrain the rest of parameters. Fits also discard the first and last three points of 
the spatial separation range. 

Simulations were performed using an Iwasaki lattice gauge action, on a lattice of dimensions 
24 3 x 6 and /3 £ [2.35, . . . ,2.55] in order to be in the range of temperatures [249, . . . , 347] MeV. The 
phase transition was estimated by the Polyakov loop behavior to be at T = 288 MeV. Except for 
one run at Q = 1 all configurations were generated at Q = 0. 

The outcome for the topological susceptibility measurements is shown in figure 1 . Our results 
are shown as black dots. They follow with a good accuracy the reference results by the Regensburg 
group [11], obtained with direct zero eigenmode counting. In order to cross-check the methodology 
we also accumulated configurations without fixing topology, by using the same HMC algorithm. 
The parameters for this run were chosen so that it corresponds to the same temperature as the 
j3 = 2.50 run, T/Tc ~ 1.1. We then extracted two subsets of configurations with Q = and 2=1 
and applied the correlator methodology to extract topological susceptibility. The final results agree 
among two sectors and moreover are in accordance with an estimate using the Atiyah-Singer's 
index theorem for all configurations (points indicated by cross Q =all, blue dot 2 = and star 
Q = 1 in the figure). This result indicates that the term in the action that forces fixed topology 
is not preventing appearance of local topological fluctuations, like pairs of caloron-anticalorons 
(so called topological molecules) which are associated to very low eigenmodes. These topological 
fluctuations are distributed in the same way as the ones generated without fixing topology. 

So far we demonstrated that our method works very well even in the finite temperature case. 
However we discovered and reported at the Lattice conference that topology fixing poses potential 
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Figure 1: (Topological susceptibility on the 24 3 x 6 lattice. Statistical errors only. Reference data points are 
diamonds and squares. The black dots are our results from measurements at Q=0. 
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Figure 2: Lowest eigenvalues at several temperatures. On abscissas the HMC trajectory number, one point 
every 100 trajectories except for RUN-II where we saved one configuration every 50. 



problems in the high temperature regime, at finer lattice spacing. The lowest eigenmodes (main 
source of the topological susceptibility signal) take longer time to thermalize than usual thermody- 
namical quantities. This is a known problem but in our case the effect appeared more strongly. The 
points at j8 = 2.50 and /3 = 2.55 were the most problematic, exhibiting some dependence on the 
initial configuration, see Figure 2. In particular the run tagged RUN-II, obtained by starting from a 
thermalized configuration at /3 = 2.50 and heated up to j8 = 2.55, was giving different results from 
a similar run starting from unit configuration. This was a severe problem, mining the reliability of 
higher temperature runs. We thus further investigated the issue by accumulating more configura- 
tions, especially for the highest temperature point. Final result is shown in Figure 2. We recall that 
all runs started from a unit configuration, high temperature, and then thermalized (except RUN-II). 
Thermalization is monitored by plaquette and Polyakov Loop as standard observables. By looking 
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at the lowest eigenvalue history, the first three panels do not show any particular problem in ther- 
malization, except the expected increasing autocorrelation time going to finer lattices. With newest 
data no anomalous behavior is detected also for /3 = 2.55. At Lattice 2011 the history was much 
shorter, less than 10k trajectories, and no jumps were present. Transition between two different 
values of the lowest eigenmode is sharp but the tunneling rate is non-zero, and no dependence on 
the starting configuration is observed. The issue was just statistical and is not expected to affect 
full QCD results. 

We thus conclude that physical quantities can be obtained from fixed topology simulations at 
finite temperature. 

3. Dynamical overlap simulations 

In this section we will discuss some of the results from simulations of two flavors of dynamical 
overlap fermions in the high temperature region. 

We accumulated 0(200-300) thermalized configurations per j8 with several masses and tem- 
peratures, see table 1. Accumulating the corresponding zero temperature configurations would 
have been a really demanding task so we do not have currently a measurement of the lattice spac- 
ing and pion mass for every j8. At j8 = 2.30 and bare quark mass am = 0.015 was 286 MeV, 360 
MeV at am = 0.025 [12]. By looking at the Dirac operator spectral density we observe that a gap 
is already opened at a temperature of around 192 MeV (/3 = 2.25) near the chiral limit, am = 0.01, 
confirmed by higher temperature simulations. This could be a signal of strong suppression of axial 
symmetry breaking, because the splitting between flavor-singlet and non-singlet correlators domi- 
nantly comes from the near-zero modes, as we demonstrate below. We do not observe this behavior 
with current data at the temperature of 177 MeV (/3 = 2.20, am > 0.025). 

We report also our first observation of the meson correlators in this range of temperatures 
(figure 4). We plot the singlet and non-singlet pseudo-scalars [t]' and n) and singlet and non- 
singlet scalars (a and 8). 

Data at jS = 2.20 do not show any degeneracy among the correlators of the lightest mesons. 
This could be an effect of the high mass (am = 0.025) but it is also possible that this temperature 
is still in the chirally broken region (see the spectral density plot). The number of configurations 
is not sufficient to precisely determine the transition temperature, so we cannot settle this question. 
The meson correlators at the highest temperature, starting from /3 = 2.25, show a clear tendency 
toward degeneracy in the chiral limit, a signal that U(\)a axial symmetry is effectively restored. 

We need to accumulate more data in order to be conclusive, but certainly we observe robust 
signals of axial symmetry restoration. Further study around the phase transition is on the way. 

4. Conclusions and perspectives 

We have shown that it is possible to perform finite temperature lattice simulations with overlap 
fermions at fixed topology. It is tested in the case of pure gauge theory where measurements of 
topological susceptibility can be compared with known results. 

Having under control the systematics of topology fixing we performed simulations of two 
flavors of overlap fermions at several temperatures in the range [170 ~ 208] MeV. We measured the 
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Figure 3: Spectral density of the overlap Dirac operator for Nf = 2. The several j3s were isolated to 
emphasize the mass dependence of the density. Zero counting of eigenvalues is intended on the left when 
line stops. 

correlators of meson operators in all the (pseudo)scalar channels looking for their degeneracy in the 
chiral limit, signal of effective restoration of both chiral and axial symmetry. We found evidence of 
this restoration, corroborated also by the spectral density analysis that exhibits a gap in the chiral 
limit at temperatures above > 192 MeV, in the current data set. 

We can currently fairly say that we have clear evidence of U(1)a effective restoration in a 
region just above the chiral phase transition in two flavors QCD. The next step is narrowing the 
region of uncertainty about the temperature when the gap starts opening. In a forthcoming paper in 
preparation a complete analysis and the newly collected data will be presented. 

This work is supported in part by the HPCI Strategic Program of Ministry of Education and 
in part by the Grant-in-Aid for Scientific Research on Innovative Areas (No. 2004: 20105001, 
20105003, 20105005, 21674002, 21684013). 
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0.144 


171 


0.05 


2.20 
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Table 1: Parameters of some of the collected data for two flavors of dynamical overlap fermions, lattice 
dimensions 16 3 x 8. The reported lattice spacings are the extrapolations to the chiral limit. 
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Figure 4: Meson correlators at several temperatures and masses. 
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